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We argue that ihe freezing transition scenario, previously explored in the statistical mechanics of 1 //—noise 
random energy models, also determines the value distribution of the maximum of the modulus of the character- 
istic polynomials of large N x N random unitary (CUE) matrices. We postulate that our results extend to the 
extreme values taken by the Riemann zeta-function (^{s) over sections of the critical line s = 1/2 + it of con- 
stant length and present the results of numerical computations in support. Our main purpose is to draw attention 
to possible connections between the statistical mechanics of random energy landscapes, random matrix theory, 
and the theory of the Riemann zeta function. 
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The Riemann zeta-function 

n=l p ^ ^ ^ 

encodes the distribution of the primes p in the positions of 
its non-trivial zeros The Riemann Hypothesis, which is 
of central importance in mathematics, places these zeros on 
the critical line s = 1/2 + it, t e M. There has in recent 
years been considerable interest, stemming from speculations 
about a spectral interpretation of the Riemann zeros, in pos- 
sible connections between the theory of the zeta-function and 
quantum mechanics (see, e.g. [2|). Our focus here is on sug- 
gesting a new link with physics via the statistical mechanics 
of disordered landscapes. 

Some of the most significant questions in the theory of the 
zeta function concern the distribution of values it takes on the 
critical line, where it behaves like a quasi-random function of 
t (essentially because the phases of the terms in the product 
in ([TJ contribute as if they were random) . It was proved by 
Selberg, for example, that, as t oo, log|C(l/2 + it) \ be- 
haves like a Gaussian random variable with mean zero and 
variance ^ log log see e.g. Q. In other words, the typ- 
ical size of log I C( 1/2 + it)\ is of the order of -y/log log t. 
As regards the extreme values taken by the zeta function over 
long ranges, the Lindelof Hypothesis asserts that |(^(l/2 + ii)| 
grows more slowly than any power of < as i — t' oo, the Rie- 
mann Hypothesis implies that it is 0(exp(const./(i))), where 
f{t) = logi/loglogi, while it is kno wn th at |C(l/2 + it)\ 
takes values at least as large as exp( ^//(i)) infinitely often 
lU pp. 354, 209]. The extreme values thus lie in the range 
between these upper and lower bounds (and are significantly 
larger than the typical values). The problem of determining 
where precisely within this range they lie has attracted con- 
siderable attention, but remains unresolved. The extreme val- 
ues in question are so rare that numerical computations have 
failed to settle the matter so far. 

One model, proposed by Montgomery, that leads to predic- 
tions for the extreme values is based on the assumption that 



the local maxima of log |C(l/2 + it)\ are statistically inde- 
pendent and that they obey the normal distribution proved for 
general values. Under such assumptions one finds that the 
typical size of the maximum value of |C(l/2 + should be 

of the order of exp (^cm \/log{t) log log(t)^ H, where cm is 
a constant. This is closer to the lower bound than the upper 
bound, implying that the extreme values are not much larger 
than the largest value known to be reached infinitely often. 

However, the values of log |C(l/2+it)| are, in fact, strongly 
correlated in a way that is significant from the point of view 
we shall explore here. Specifically, define, for a fixed t e M, 

vl'^\x)^~2\og\c(^^+z{t + x)y (2) 

Selberg's theorem implies that v}'''^ (x) behaves like a Gaus- 
sian random function of x when t —?' oo. Such a random 
process is characterized by the two-point correlation function 

C(a;i,a;2) — v/'"'' (xi)v/^-' (X2), with the bar denoting the av- 
erage over intervals [t — h/2,t + h/2], with h ^ t chosen 
so that the intervals contain, asymptotically (as t -> 00), an 
increasing number of zeros. A straightforward calculation, 
sketched in ||5l, see also 0, leads to 

^, r -21og|a;i -X2I, for < |a;i - xzl < 1 

2 log log t, for |a;i-:E2|«j^ 

(3) 

The significance of the logarithmic form of the correlations 
will become apparent when we come to make comparisons 
with corresponding problems in random matrix theory and 
statistical mechanics. 

Over the past 40 years, inspired by the pioneering work of 
Montgomery fT\ and Odlyzko |8|, considerable evidence has 
accumulated for connections between the statistical properties 
of the Riemann zeta-function and those of large random ma- 
trices. For example, correlations between the nontrivial ze- 
ros of the zeta function are believed to coincide asymptoti- 
cally with those between the eigenvalues of large random uni- 
tary or hermitian matrices lilnSj, and the value distribution of 



2 



C(l/2 + it) is believed to be related to that of the characteris- 
tic polynomials of these matrices 19r il2i . Our purpose here is 
to connect these two areas of research to a third, the statisti- 
cal mechanics of disordered logarithmically-correlated energy 
landscapes, see ifTSHTSjl and references therein. The analogy 
we develop suggests that the spin-glass-like freezing transi- 
tion that dominates the low-temperature behaviour in the sta- 
tistical mechanical problems also governs the extreme values 
taken by the characteristic polynomials of random matrices 
and the zeta function. This sheds new light on the longstand- 
ing problem of determining the maximum size of ({1/2 + it). 

To establish such a connection we consider the ensemble of 
NxN unitary matrices, chosen uniformly at random from the 
unitary group U{N) (i.e. from the Circular Unitary Ensemble, 
or CUE). We denote the eigenvalues of a given matrix l/jy by 
exp . . . , exp {i(pN) and the corresponding characteris- 

tic polynomial by 



N 

Pn{9) = det (1 - Un e-^'^) = U e-'^^"-")) 



(4) 



n=l 



It is instructive to compare v}'^\x) from ^ with vl;j^\9) = 
-21og |pAr(6')|. It was proved in |9| that vl:J^\9) satisfies 
a central limit theorem: the distribution of the values of 
log |pAf(0)| tends to a normal with mean zero and variance 
21ogA^ as iV — > oo. Identifying the mean density of the 
eigenvalues, N/2tt, with the mean density of the Riemann ze- 
ros near to height t, log renders the agreement with 
Selberg's theorem for log |C(l/2 + it)\ complete. Importantly 
for us here, V^^^ has the following representation 1 10|: 



-me (AT) 



comp. conj. 



(5) 



with Vnwi^^ = Tr(C/^")). According to Q3, as Un 

varies in the CUE the coefficients v'^^ for any fixed n 
tend in the limit iV — > 00 to independent, identically dis- 
tributed complex gaussian variables with zero mean and 
unit variance. Denoting averages over the unitary group 
U {N) with the angular brackets, a simple calculation shows 

that (l4^^(6'i)l4^'(6'2)^ tends in the limit N 00 

to — 21og (2| sin ^{Oi — ^2)!), and so exhibits precisely the 
same logarithmic behaviour at small distances as in ([3]). For 
large but finite N, the logarithmic divergence can be shown 
to saturate at |0i — (?2| ^ N^^, so after associating N ^ 
log ^ the correspondence between the small-scale behaviour 

of vj^"^ (9) and v/'''' (x) becomes complete. This is significant 
from the point of view we seek to develop. 

Our primary goal is to determine the distribution of the 
maximum value of jpTvl^*)! over < < 27r when the ma- 
trix Un ranges over U {N). Our second goal is then to use the 
random-matrix results to motivate predictions for the extreme 
values of the Riemann zeta function. The first steps in this 
direction were taken in [4], where the tail of the distribution 



of the maximum values of \pi\f{9) \ in < 9 < 2-k was found. 
This tail determines the typical size of the maximum values 
expected when C/^v is sampled independently a large (expo- 
nentially in N) number of times from within U{N). The re- 
sults are consistent with the Montgomery model and lead to a 
prediction for the value of the constant cm there. 

The focus here will differ from that of |4| in the following 
ways. We shall be concerned with the distribution of max- 
imum values of the characteristic polynomials of single ma- 
trices, rather than with large numbers of matrices, and will 
obtain the full value distribution of the maxima in the limit as 
N 00, rather than concentrating on the tail that is relevant 
when maximizing over many matrices. This leads to a model 
for the distribution of maximum values of |C(l/2 + it)\ over 
the intervals T < t < T + 2tt, rather than < t < T as 
T 00. Furthermore, it makes the problem of numerical 
computation of the distribution in question significantly eas- 
ier, because one is finding the maximum of ^ logT rather 
than ^ TlogT numbers (the values of the local maxima). 

We start by noting that the maximum value ofp]y{9) can be 
characterized in terms of 



ZnW) 



N 
2tt 



27r 



\pM{9)\'^d9^ 



N 
2ti 



27r 



e-^^"(^) d9, 



(6) 

where Vn{G) = -21og |pAr(6l)| and /3 > 0. SpecificaUy, if 
= -r' log then 



lim TiR) = min Vn{.9) 

P^oo ee(0,27r) 



-2 max log |pjv(6')|. 

0e(O,27r) 



(7) 

The key observation is that (j6]) takes the form of the partition 
function for a system with energy Vn{9) and inverse temper- 
ature (3. J- {(3) may then be associated with the corresponding 
free energy. Recalling that the values of Vn{9) are asymptot- 
ically gaussian distributed and logarithmically correlated, it is 
natural to draw comparisons with a class of problems that has 
attracted a good deal of attention recently in the area of disor- 
dered systems, namely the statistical mechanics of a particle 
equilibrated in a random landscape with logarithmic correla- 
tions. Two-dimensional systems of that sort appear in many 
contexts, for example in the problem of Dirac fermions in a 
random magnetic field ifTTI . and one-dimensional analogues 
were considered recently in lfT4l [TS l. Importantly for us, in 
the one-dimensional case the energy landscape is given by a 
regularization of the random Fourier series featuring in Q. In 
the statistical mechanical problem there has been a particular 
focus on the so-called freezing transition [13.. 17] which dom- 
inates the low temperature limit and determines the extreme 
value statistics. The latter appears to be manifestly different 
from the case of uncorrelated variables. We shall argue that a 
similar freezing transition determines the extreme value statis- 
tics of the characteristic polynomials and hence, conjecturally, 
of |C(l/2+it) |, and that therefore the logarithmic correlations 
exhibited by V^'^'' (9) and v/''-' {x) play an important role. 

The first step is to consider the positive integer moments 
(^Z^{(3)) , k = 1,2, . . .. Using standard methods of ran- 
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dom matrix theory (RMT) it is easy to show that 



---Jo Df{Q)fA^^ 



where D 



(/?) = detfAi^f')^ ' is 



Toeplitz matrix 



M;^{- = / e'-^^'-J^ TT [2 -2 cos I 





the determinant of a 



0,)f'l- (9) 



27r 



When — > oo the asymptotics of such Toeplitz determinants 
is well known [TSl to be given by 



D 



N 



G(l + 2/3) 



k 

n 



-2,92 



(10) 



where G(x) is the so-called Barnes function. Substituting ( 10 1 
back to (j8]l we see that the resulting expression is the standard 
Dyson-Morris version of the Selberg integral [19J , convergent 
for k < -02, and divergent for larger k. As we have A: > 1 
the procedure makes sense only for < 1. Introducing 



G(l+2/3)r(l-^i2)' 

(Z^(/3)) = Z^T{1 



we find 



kl3^), k<p- 



(11) 



The expression ( [TT] ) for the moments has exactly the same 
form as that for the partition function of the landscape of the 
circular-logarithmic model (a periodic version of 1 // noise) 
lfT4l . but with a different value for the characteristic scale Z^,. 
This means we can simply translate, mutatis mutandis, the 
results of |14 | to the values of characteristic polynomial sam- 
pled along the full circle 6* e (0, 27r]. In particular, we con- 
clude that the maximum value of the modulus of a CUE char- 
acteristic polynomial pn{0) in an interval 9 g [0, 2tt] can be 
written in the limit iV — > 00 as 



2 max log|pAr(0)| 



(12) 



where ajv = — 21ogA^ + clog log iV + o(l), with, conjec- 
turally, c = |, ^at = 1 + 0(1/ log A^), and a: is a random 
variable taking values distributed with probability density 

pix) = 2e-ifo(2e"/'), (13) 

where K,y{z) denotes the modified Bessel function. Two con- 
sequences are particularly noteworthy: (i) the tail of the proba- 
bility density of the random variable x is given asymptotically 
when X — > —00 by p{x) ~ \x\e^, (ii) the value of the con- 
stant c = 3/2. These two features are believed to be universal 
characteristics of the extreme value statistics of the class of 
logarithmically correlated random variables |13 |, and distin- 
guish those from short-range correlated random variables, for 
which c — 1/2 and p{x) ~ (Gumbel distribution) |20|. 
This new class is believed to include, in particular, the 2D 



Gaussian free field II2TI . branching random walks ll22ll . poly- 
mers on disordered trees ||23l , 1// noise lfT4l[T5l and models 
appearing in turbulence and financial mathematics ll24l . 

In order to illustrate the extreme value predictions in the 
context of the Riemann zeta function we now summarize 
the results of preliminary numerical computations. These 
involved numerically evaluating C(l/2 + it) over ranges of 
length 27r, at various heights T, and finding the maximum 
value Cmax(27r; T) in each range. We used the amortized- 
complexity algorithm [25 1, which is suitable for computing 
C(l/2 + it) at many points. 

The first test concerns the value of the constant c in ( [T2] i. 
We expect the logarithmic correlations to lead to c = |, rather 
than c = 5> would be the case if the zeta correlations 
were short-range. The mean of Cmax(27r; T) suggested by the 
model in ([12]) is (5 = e'^ N/(logN)^ , with c = 1/2 or 3/2, 
and 7 — 0.57721 . . . , where we set N to be the nearest in- 
teger to logT. At each height a sample that spans « 10^ ze- 
ros is used yielding « 10^/A^ sample points (since there are 
roughly N zeros in each range of length 2tt). In view of the 
table below, one may conclude that c = 3/2 fits the data con- 
siderably better than c = 1/2, thus supporting the logarithmic 
correlations model. 



TABLE I: Ratio of data mean S to model mean S with c = 3/2 and c = 1/2. 



T 


TV 


(~5/5) 

V /c=3/2 


V /c=l/2 


1022 


51 


1.001343 


0.504993 


1019 


44 


0.992672 


0.510293 


1015 


35 


0.976830 


0.518057 


3.6 X lO'' 


17 


0.930533 


0.552856 



Testing the distribution p{x) is more difficult, because 
the data converge extremely slowly at that scale. The re- 
sults of our initial experiments are summarized in Figure [T] 
Specifically, we consider — 21og |Cinax(27r;r)| + 21ogA^ — 
I log log based on a set of approximately 2.5 x 10* ze- 



ros near T 



10 



28 



The data are normalized so that 



-2 log |Cmax(27r; T) I + 2 log iV - I log log N has empirical 



/ x'^p{x) dx 



3.28986813. 



The overall 

agreement is supportive of ( [T3] l, especially in the important 
tail when .t — > — 00 and in view of the fact that lower order 
arithmetical terms |9 , 1 1 1 have not been incorporated, but can- 
not be said to be a conclusive verification at this stage. The 
behaviour in the tail is significant because if it were to per- 
sist into the large deviation regime it would suggest that the 
Montgomery heuristic significantly underestimates the max- 
imum values achieved by |C(l/2 + it)\, although this seems 
unlikely. 

Finally, we put the specific idea of the freezing to a numer- 
ical test by considering the analogue of the partition function 
Zisi{li) defined by (j6| which is clearly 



Nt 
2^ 



T+2ir 



|C(l/2 + Z2/)p%, 



(14) 
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FIG. 1: Numerical computation (solid red line) compared to theoret- 
ical prediction i 13 i (dashed black line) forp(3;). 



data 
freezing 



FIG. 2: Numerical computation (red dots) compared to the theoret- 
ical prediction (dashed black line) for — /(/3), suggesting freezing 
beyond /? = 1 



where Nt = log and calculating the ensuing free energy 
/(/3) = —{/3 log (iVT))~^log zp by averaging over 10^ values 
of T near T — 10^®. If freezing is operative, -/(/3) is expected 
to be equal to /3 + 1//3 for /3 < /3c = 1 and remain frozen to 
= 2 for all ^ > 1 ifllUBjl . The results shown in Figure 
2 (which are renormalized by the subtraction of lower order 
arithmetic terms |5|) again are consistent with this prediction. 

To conclude, we have put forward speculations on the ex- 
istence and implications of the freezing phenomena seen in 
the statistical mechanics of random energy landscape models 
in the context of the extreme value distribution of the charac- 
teristic polynomials of large random matrices and, further, of 
the Riemann zeta-function on the critical line. We believe this 
sheds interesting new light on these problems. In a forthcom- 
ing paper [5] we present the details of our calculations and 
discuss further the broader picture of freezing phenomena and 



their manifestation in the present context. 
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